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Abstract 

We show that the complex "PT-symmetric periodic potential V{x) = — (i^sin2x-|- 
A^)■^, where is real and is a positive integer, is quasi-exactly solvable. For odd 
values of A > 3, it may lead to exceptional points depending upon the strength of the 
coupling parameter ^. The corresponding Schrodinger equation is also shown to go 
over to the Mathieu equation asymptotically. The limiting value of the exceptional 
points derived in our scheme is consistent with known branch-point singularities of 
Mathieu equation. 
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In recent times, considerable attention has been paid to the development (see, for exam- 
ple, the comprehensive review [Ij) of complex non-Hermitian systems admitting space-time 
reflection [VT) symmetry, the primary reason being rooted to the conjecture that the eigen- 
values of the governing Hamiltonians support [2] a real bound-state spectrum. However, the 
package of PT-symmetric theories is such that should a spontaneous breaking take place, 
a loss of the reality of the energies is incurred and complex eigenvalues appear in conjugate 
pairs. In such a situation, the accompanying eigenfunctions cease to belong to the VT 
operator. The critical points where some pairs of real eigenvalues transit to the complex 
plane turn out to be exceptional points, i.e., points of the complex parameter space where 
two eigenvalues coalesce [3l Occurrences of the latter have been widely studied in the 
literature (see, for example, [5], El [Tl E] ) . 

This note investigates a new class of quasi-exactly solvable (QES) PT-symmetric vr- 
periodic potentials defined by 



with ^ G M and an integer, which has applications in one-dimensional crystal problems [9] . 
We shall be interested here in vr-periodic or vr-anti-periodic solutions of the corresponding 
Schrodinger equation, which are counterparts of the band edge wavefunctions (often called 
eigenfunctions) of conventional theories for real periodic potentials (see, for example, ref. 
[To]). As will be demonstrated below, for odd values of (> 3), V{x) leads to complex- 
conjugate square-root branch-point singularities, implying the crossing of the energy levels 
to the complex plane and giving rise to exceptional points. Here the eigenfunctions turn 
out to be PT-symmetric and vr-periodic in nature. In contrast, for even A^, such critical 
points are absent: there are only pairs of complex-conjugate eigenvalues with corresponding 
eigenfunctions that are vr-anti-periodic and fail to be PT-symmetric. 
It should be noted that the 'hyperbolic' analogue of ([1]), namely pT] 



^It should be noted that the critical points considered here are only special cases of exceptional points 
because (i) we restrict the range of ^ in ([1]) to real values, while we might consider the whole complex 
plane, and (ii) exceptional points do not necessarily occur as branch-point singularities of a single analytic 
function, but might also result from the equality of two analytic functions. 





(Csinh2?/ -iMf, 



(2) 
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where C G M and M an integer, is also PT-symmetric but always has energies which are real 
for all integer values of M and as such no critical point is encountered. A similar conclusion 
holds for the cosh- version of ([2]) [12] (see also [13]) whenever M is an even integer because 
the QES eigenvalues occur in complex-conjugate pairs, but it does not hold for M an odd 
integer since then the QES eigenvalues are real for | C, \ smaller than or equal to some critical 
value. 

Calculations carried out by us show that for the QES eigenvalues of ([1]), critical points 
emerge when the quantity A^|^| keeps to a value in the vicinity of 1.47 for G Nodd- This 
prompts us to enquire into a large behaviour of the Schrodinger equation having V{x) 
as its potential. In units wherein h = 2m = 1, the Schrodinger equation reads 



dx 
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{i^sm2x + Nyijj = Eijj (3) 



and it is a simple exercise to convince oneself that in the large A^ limit, equation ([3]) goes 
to the form 

- ----2ig sin 2xip = E^), (4) 
dx^ 

where g = G M and E = E + N"^. In writing (jlj), we have assumed g to be finite under 
the double limits N oo and ^ ^ 0. 

Interestingly, equation (jlj) is the Mathieu equation [14J, in which branch-point singu- 
larities are known to occur as exemplified in the works of Blanch and Clemm [15], Hunter 
and Guerrieri pjj, and others [T7l [TSl [T9]. A variety of numerical methods has attempted 
to solve Mathieu equation with an imaginary characteristic parameter as is the case in 
(jlj). Tabular values from an old paper of Mulholland and Goldstein pTT] reflect^ that the 
eigenvalues E^ would be both real and conjugate complex: the transition occurring at 
g = gc = 1.4687 for solutions having a period vr. It is of utmost relevance that gc is consis- 
tent with our calculated transitional value of A^|^| ~ 1.47 (where A^ is odd corresponding to 
TT-periodic eigenfunctions) beyond which the energy degeneracy is lifted, the critical points 
effecting a phase transition. 



^ after making a change of the coordinate x ^ t — j without affecting the underlying periodicity. 
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In a different context, Bender et al [20] have made a systematic WKB analysis of the 
potentials of form V{y) = i sin^""*"^ (y) (n = 0, 1, 2, . . . ) to justify that these have real 
bands and gaps, the wavefunctions turning out to be always periodic at the edges of each 
band. The presence of the coupling g in the Mathieu equation (jlj), however, makes it of a 
more general character and facilitates the study of critical points. 

We now turn to and proceed to interpret it as a QES system. Typically QES 
problems may be generated by means of a universal enveloping algebra U(sl(2)) [2TI [22] . 
Here we observe that in terms of the sl(2) generators 

J+ = z^-i- - 2j2;, Jo = z-^ - j, J_ = (5) 
az dz az 

the Hamiltonian H can be mapped to the following spectral problem 

ifg = 4 J2 + 2^J+ + 2^J_ -N^ + ^^ (6) 

for j = {N — l)/2, where in principle both even and odd values of N are allowed. To arrive 
at the representation ([6]) we utilized the connections 

H, = [f^izT'Hfiiz], 

M^) = .(^-^)/^exp f.-i 
4 V z 



(7) 



with z = exp(2ix), as is required to have a gauge-transformed form. 

The eigenvalues and eigenfunctions of H can now be derived from those of ([2]) by means 
of an anti-isospectral transformation y ^ ix, ( ^ — i^, M ^ N [231 El]- We can check 
that the gauge function and the gauged Hamiltonian written down in [11] transform to fi{z) 
and —Hg as they should. Thus equations ([6]) and ([7]) are in the true spirit of quasi-exact 
solvability. 

Our results for some of the eigenvalues E and eigenfunctions iplx) of H are furnished 
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below (cr, r = ±): 

mN=l: E = -l + ^^ oc es^^'-^^x^ 

• N = 2: E„ = -3 + 2ia^ + ^'^ , ip^{x) oc e^^'''''^''{cosx - a sinx), 

• N = 3: Eo = -5 + ^'^, ^/'o(x) oc es^'^^'^^x 2x, 



2isin2x+ I (^1 + (Tv^I - 4^^ 
AT = 4: E,,, = -ll-2iae + e'-4rv/rTT^e^, 



i/jfj^rix) OC e2€'^™^^'(cosx + (Tsinx) isin2a; + | ( 1 + Ta/I + icr^ — ^ 



The expressions for the energy and the associated eigenf unctions for higher values of N 
are more cumbersome: of the five solutions for N = 5, one set of energy eigenvalues 



is obtained by solving a quadratic equation and the corresponding eigenf unctions are 

1 



(9) 



Mx) oc e5«'^''^2x^Qg2a; 
while the other set corresponds to the cubic equation 



2isin2a; + - (3 + (7^9-4^2 



(10) 



£^ + 20£^ + 64(1 + ^2)^ + 768^2 ^ 
where S = ^"^ — 25 — E, with eigenf unctions 

^p{x) oc e^^''''^''(sin^2x + 6sin2x + c). 
In (fT2l) . the quantities b and c are given by the relationships 



:ii^ 



(12) 



(13) 

The three solutions of £ from (fTTil . when substituted in (fT3l) . yield the corresponding values 
of b and c which, in turn, provide the expressions for the eigenfunctions from ( |T2l) . 

We thus notice that for even A^, VT symmetry is broken: only pairs of complex- 
conjugate eigenvalues exist, the corresponding eigenfunctions being vr-anti-periodic. On 



the other hand, for odd A^, there is at least one real eigenvalue irrespective of whether VT 
symmetry is broken or not, while the remaining ones are either real or in complex-conjugate 
pairs depending upon a critical value of ^. However, the corresponding eigenfunctions are 
TT-periodic but their being PT-symmetric or not depends upon the nature of eigenvalues. 

While it is trivial to see from ([8]) that for = 3 there are critical points at |^| = 
0.5, determining those for the case = 5 requires solving the cubic equation ffTTl) . A 
little algebra shows that the condition for the existence of one real and a pair of complex- 
conjugate roots relies upon the positivity of A that is defined by 

A = 16^^ - 4^^ + 103^2 - 9. (14) 

Solving for A = gives the roots ^ = 0.0876 and 0.0812 ± i(2.5331), from which we 
immediately identify critical points in the neighbourhood of |^| = 0.296. The other branch- 
point singularities at |^| = 1.5, where the eigenvalues (Q coalesce, are situated further away 
from those coming from the cubic equation ffTTj) and for the purpose of comparison with 
Mathieu equation do not interest us. 

The above procedure can be continued further. For A^ = 7, we are led to a cubic 
equation 

+ 56£^ + 16(49 + 4^^)^ + 758(3 + 2^^) = q (15) 

and a quart ic one 

£^ + 56£^ + 16(49 + 10^2)^2 ^ ^q^q _^ 17^^)£ + 2304^^(24 + = 0, (16) 

where £ = ^"^ — 49 — E. A search for their zeros tells us that there are critical points 
coming from the latter and situated at |,^| = 0.2107 further which the eigenvalues develop 
square-root singularities. The corresponding eigenfunctions are too complicated to be given 
here. 

Putting everything together, it is significant to note the convergence of the quantity 
A^|.^| to the value of 1.47. Indeed our findings are 

A^l^l = 1.5 for A^ = 3 

= 1.4797 for A^ = 5 (17) 
= 1.4749 for A^ = 7 
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which are comfortably close to Mulholland- Goldstein estimate of gc = 1.4687. The trend 
prevails for higher odd integers as well, as our numerical results have indicated. We 
therefore conclude that the QES model given by ([T]) is well described by an imaginary sine 
potential in the large N limit. The presence of the imaginary coupling constant i^ makes 
the role of the critical points very specific as comparison with the numerical solutions of 
the Mathieu equation has shown. 

Summarizing, we have proposed a QES PT-symmetric periodic potential which pos- 
sesses real and/or complex-conjugate eigenvalues. The latter may assume square- root 
branch-point singularities signalling the appearance of exceptional points. It is also es- 
tablished that our model goes over to the Mathieu equation asymptotically and that, even 
for rather moderate odd values of A^, there is a remarkable coincidence of the limiting 
value of the lowest exceptional points found in our scheme with the estimates of the lowest 
branch-point singularities of Mathieu equation [T71 118] . 
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